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ABSTRACT 


The principal purpose of this study is to investigate the feasibility 
of using the Scanning Celestial Attitude Determination System (SCADS) 
during Earth Resources Technology Satellite (ERTS) missions to compute 
an accurate spacecraft attitude by use of stellar measurements. The 
spacecraft is that designed by G. E. Corporation, and is local-vertical- 
stablized. 

A heuristic discussion of the SCADS concept is first given. Two concepts 
are introduced: a passive system which contains no moving parts, and an 

active system in which the reticle is caused to rotate about the sensor *s 
axis. Later it is shown that the aperture diameter may be one inch for 
the passive system, and two inches for the active system. 

A quite complete development of the equations of attitude motions is 
then given. These equations are used to generate the true attitude which 
in turn is used to compute the transit times of detectable stars and to 
determine the errors associated with the SCADS attitude. 

The stray light problem which originates because of direct sunlight, 
earth-reflected sunlight , and spacecraft-reflected sunlight is then 
/ considered. It is shown that a location aboard the spacecraft's lower 

sensor ring does exist which minimizes the stray light problem and also 
provides a favorable location for the SCADS sensor. 

A more complete discussion of the analytical foundation of SCADS concept 
and its use for the geometries particular to this study is then given. 
Also given are salient design parameters for the passive and active 
systems . 

An error analysis indicates the errors of the passive and active systems 
(two head sensor) to be on the order of 0.04^ and 0.02°, respectively. 
These errors are associated with a worst-case attitude control and a 
worst-case stellar background. 

Finally, suggestions are given which might reduce these attitude errors. 
The suggestions involve changes in the slit-pattern and method of data 
reduction . 
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INTRODUCTION 

Consider a nearly local-vertical- stabilized satellite orbiting the 
earth. Let denote the local-vertical^ system of coordinates which is 
associated with the instantaneous satellite position. Also, let be a 
system of coordinates fixed in some portion of the satellite (S^ and 
will be defined more precisely in the next sections) . The orientation of 
with respect to may be defined by three angles ^ , 0 , cp which will 

be called yaw, pitch, and roll, respectively. The general problem is to 
determine the satellite's yaw, pitch, and roll at any given time by utilizing 
a particular set of onboard stellar measurements. 


Suppose an optical system is constructed so as to provide a set of 

/V. /V 

viewing directions which lie on a curve d(T]) , T|^ < T| < T\^. Where d(T|) is 
a unit vector whose components with respect to are known if T\ were given. 

Since S^ is changing orientation with respect to an inertial system, 

A 

the curve of viewing directions d(T|) will intermittently sweep across a 
star. This event is called a transit and the instant of its occurrence is 
called a transit time . 


At any transit time, which is produced by a known star, the spacecraft’s 
attitude is constrained. This constraint is generated as follows: 

(1) The direction to the star resolved in is known except for 
a single parameter, T], 

(2) This direction may be resolved in by introducing yaw, pitch, 
and roll at the measured time. Thus, the direction resolved in 

is known except for four parameters: T|, \(f(t.), 0(t.), cp(t.), 

where t^ is the i^^ measured transit time. 


It is assumed that the satellite’s position is a known function of time; 
hence, the orientation of with respect to a preferred inertial system 
is known. 

2 

The problem of identifying the transited star is not trival. However, it 
will not be considered here. A discussion of the problem and its solution 
for two particular geometries may be found in References 1 and 2. 
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(3) Since the star is identified, its direction with respect to 

is readily computable with the aid of the star catalog and the 
satellite’s known position, 

(4) Equating these two directions to the star yields two independent 

equations in the unknowns T|, 0(t^), cp(tp . 

(5) The parasitic unknown, T), may be eliminated from these two 
equations to obtain one equation in the unknowns \|f(t,), 0(t.), 
and cp(t^) . 


The equation generated at Step (5) will be called the constraint 
equation implied by the measured transit time. The ease with which T| may 
be eliminated from the two equations obtained at Step (4) is dependent 
upon the functional form of d(T|). The elimination is particularly simple 
if the directions d(T|) lie on a right circular cone. Fortunately, the 
instrumentation of forming d(T|) to lie on such a cone is simple. Even 
more simplification is obtained if d(T|) lies in a plane. This type of 
configuration is assumed in this study, and is pictured schematically 
in Figure 1. Because of the instrumentation, we call the optical element 
which generates the viewing directions, d(Tl) , a slit; and the portion of 
a plane which contains d(7]) a slit-plane. 

Two problems, which are given considerable attention later in the study, 
may be noted. The first may be called the transit-rate problem; the second, 
the modeling problem. 

Transit Rate 

Our discussion thus far defined a slit-plane to have a fixed orienta- 
tion with respect to S^. Such a system will be called a passive system. 

At any instant of time, a star may or may not lie in this plane. If a 
star is not detected, no information is obtained. Transits will be gathered, 
however, since changes orientation with respect to an inertial system. 

But, this change in orientation is rather slow (approximately one revolution/ 
orbital period) ; so, too few transits tend to be gathered. Some methods of 
overcoming this difficulty are listed in Table I. Also listed is the 
qualitative effect of each change. 
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0 



Figure 1: Schematic of Radial Slit Instrumentation. 
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TABLE I 

METHODS OF INCREASING DATA RATE 


Parameter Change 

Qualitative Effect 

1. Increase field of view 

Adds more stars, more background noise, 
increases optical inaccuracies or adds 
to optical costs, adds more weight. 

2. Increase number of slits 

Adds more transits/star, more background 
noise, greater fabrication difficulties, 
complicates star identification , 
increases possibility of near simultan- 
eous transits (since only one photo- 
detector is used) . 

3. Increase number of sensor 
heads 

Makes for better stellar geometry, adds 
candidate stars , increases cost , s ize , 
weight , power requirements ; but , reduces 
size, etc.,^ of each individual head. 

4. Sense dimmer stars 

Adds candidate stars, complicates star 
identification, increases optical 
aperture, puts more load on optics and 
electronics in terms of accuracy and 
sensitivity, increases probability of 
false transits. 

5. Rotate the reticle 

about the optical axis 

Adds more transits/star, requires an 
angle encoder or its equivalent to 
sense orientation of reticle with respect 
to Sy , increases optical aperture, 
increases cost, decreases reliability. 
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Parameter changes 2 and 3 can be used very effectively for the 
passive system. Hence, the recommended passive system consists of two 
heads, in which each sensor head possesses 5 slits (Figures 6 and 16). 

Parameter change 5 yields a configuration which will be called an ”active‘* 
system. This term is used because the reticle is made to rotate within the 
instrument. In this case, the viewing directions are specified as 
d(T], t) or d(T|, 0), where t is time and 0 is an angle encoder reading. 

Again, d is known except for a single parameter T] so the previous 
discussion still holds for the active system. 

The advantages and disadvantages of a passive system with respect 
to an active system for a nearly local-vertical-stabilized satellite 
are given in Table II. 

A major question, which is answered in this study, is the accuracy 
with which the passive and active systems can determine the spacecraft *s 
attitude. Another important question -- What is the optimum spin rate if 
an active system is to be used? — is also dealt with. 

Model 

The second major concern is the attitude model which dictates the 
form of the yaw, pitch, and roll motion which is to be used in reducing 
the data. The need for such a model arises because the basic constraint 
simply yields one equation in ^(t^, 0(t^), cp(t^) which are yaw, pitch, 
and roll at the measured transit time. Since only one transit time is, 
in general, gathered at t^ and perhaps no transit times are gathered at 
a time at which an attitude estimate is desired; a time-dependent 
characterization of yaw, pitch, and roll must be given beforehand. This 
attitude model may then be used to couple the constraint equations. So 
as not to yield an unwieldy or ill-conditioned numerical problem, the 
functional form of the attitude model must contain only a few unknown 
parameters . 
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TABLE II 

ADVANTAGES AND DISADVANTAGES OF A PASSIVE SYSTEM 
WITH RESPECT TO AN ACTIVE SYSTEM 


1. Size . The optical aperture tends to 
be smaller since the slits scan 
slowly with respect to the stellar 
background. Also, no drive or read- 
out mechanism of the driven reticle 
is required. 


2 . Power . 

3 . Weight . 

4 . Reliability . 


5 . Cost . 

6 . Onboard Storage . Since fewer 
transits are gathered per unit 
time, less total onboard storage 
is required. 

7 . More Nearly Uniform Stellar Distri - 
bution . Since the scan is slow, 
dim stars can be detected. As the 
spin rate of the active system is 
increased, only the brighter stars 
are detected. More transits may 

be gathered per orbital period, 
but the rate tends to be less and 
less uniform as the spin rate is 
increased , 


1. Model . Since the star transits 
are gathered more slowly, the 
transits used to compute the 
attitude are necessarily spread 
over a longer time interval. 
Thus, a more complete model of 
the attitude motion may be 
required. 

Analysis . More difficult. 

3 . Data Reduction . 

4. Attitude Error . Greater RMS 
error over all cases in one 
orbital period . 


Passive Advantages 


Passive Disadvantages 
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A second model is necessary in this study, this is the "real-world*' 
model. This study uses the attitude dynamics as governed by the ERTS 
spacecraft designed by General Electric Corporation. A rather complete 
development, including all major effects, of the equations of motion is 
given. This resulting ninth order system of differential equations is 
solved numerically to obtain the real world attitude. It is felt that 
the resulting real world attitude model is a faithful simulation of 
the physical problem, except for the attitude control law. That is, 
the control system has a horizon sensor and a rate gyro. Outputs from 
these sensors are processed and commands sent to three flywheels and 
three jets. The main attitude control is via the flywheels. The jets 
are used primarily during initial attitude acquisition. In this study, 
some liberty is taken in that it is assumed that yaw, pitch, and roll 
are somehow measured (with errors) . These measured values are then used 
to command the fl}rwheels via a control function which is a linear function 
of the measured yaw, pitch, and roll and their measured rates. 

The spacecraft orbits and attitude characteristics will now be 
briefly defined. Further explanation is given in the text. 

Orbit 

One particular class of the low altitude, sun synchronous orbits 
has been assumed in this study. The orbital parameters of this class 
are given in Table III. The qualitative effect upon this study of the 
choice of orbit is as follows: 

(1) As the orbital altitude is lowered, a greater portion of 
the stellar field is obscured by the earth; and earth 
reflected sunlight becomes a greater problem. In addition, 
the time derivative of the true anomaly is greater which, 
considered by itself, implies a greater number of stellar 
transits per second for the passive system. One probably 
cannot consider a passive system for a synchronous (24 hr.) 
orbit . 
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(2) Since the orbit is sun synchronous, the sun*s direction, as 
viewed from the spacecraft, is constrained. This constraint 
can be used to avoid direct and scattered sunlight entering 
the sensors. 


TABLE III 

ORBITAL PARAMETERS 


Parameter 

Value 

a 

4563 nautical miles 

e 

0 

i 

99.08° 

Q 

152.5° + Right Asc. of sun (9:30 a.m.) 

n 

l°/day 

Period 

6196.015 sec 

Height 

492.35 nautical miles 

True anomaly rate, v 

0.0581°/sec 


Attitude Parameters 

The spacecraft is that designed by General Electric Corporation for 
the ERTS missions. It is attitude stabilized by jets and flywheels; 
however, in this study only the flywheels are used for stabilization. 

A dynamic equivalent of the spacecraft *s configuration is given in 
Figure IB.* The spacecraft is composed of five rigid sections: 

(1) main body, 

(2) solar array, and 

(3) three flywheels. 


^Figure 1 of Appendix B. 
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The solar array rotates about a fixed axis. The angle of rotation is 
chosen so as to maximize the received solar energy. Important spacecraft 
attitude parameters are given in Table IV. 

TABLE IV 

ATTITUDE PARAMETERS 



External Torque 
-4 

Maximum magnitude = 1.15 x 10 ft lbs 
Minimum magnitude =4.71x 10 ^ ft lbs 
(See Appendix D) 

Attitude 



Yaw 5 ^ 

Pitch, 0 

Roll, cp 

max. (deg x 100) 

85.2 

38.6 

42.8 

min. (deg x 100) 

-144.7 

-75.8 

-76. 1 

(See Figure 

21) 
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SYMBOLS 


principal moments of inertia of spacecraft without 
solar array 




a, e, Q, 0 ), i 


D. 

J 

H 



9 

9 

j A9 , Acp 


9, cp 


V 

T1 


$(M, S) 


average value of principal moments of total spacecraft 

orbital elements of spacecraft’s trajectory 

angular momentum of all three flywheels resolved in the 
direction i^ , j = 1; j = 2; k^, j = 3 

angular momentum with respect to an inertial system of 
the spacecraft 

components of a star in the coordinate system S^ 
mass of solar array 

components of external torque in the system S^ 

angle which defines orientation of solar array with 
respect to the spacecraft 

^ t h 

angle between yaw axis, , and projection of the i 

sensor optical axis in the yaw-roll plane (i^, j^) 

angle defining orientation of slit j with respect to its 
optical axis 

true yaw 

true pitch 

true roll 

error in the measurement of yaw, pitch, and roll by the 
attitude control system 

estimates of yaw, pitch, and roll obtained by use of the 
star sensor system 

true anomaly of spacecraft 

slit-curve is written as a function of this parameter 
t h 

angle between i sensor optical axis and pitch axis 

inertia tensor of the set of masses, M, computed with 
respect to the coordinate system, S 
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(JD 


ij 


angular velocity of the coordinate system S. with respect 
to the coordinate system S^ ^ 


Coordinate systems : 


10 ij 


celestial system 
local geocentric system 
spacecraft fixed system 

system associated with i^^ head and slit 


In general: 
A* 

U 

<x> 

U 


is the transpose of a matrix A 
implies U a unit vector 

is the average value of x(t) in a given time interval 
equals 

implies U a vector with three components 
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SYSTEMS OF COORDINATES 


To aid in later discussion, a set of fundamental coordinate systems, 

S., will now be introduced. A triad of orthonorraal right-handed vectors 
1 

A yv ^ 

associated with each system, S^, will be denoted by (i^, j^, k^) . Also, 
^i+1 ^i such that one of the following relations is 

satisfied : 


^i+1 ^i^ 


j . . 1 = j . > k . . = k . 
‘^i+l 1+1 1 


- the celestial system, i^ in the direction of the First Point of 
Aries, k^ in the direction of the north pole. This system is needed 
because stars have their directions catalogued with respect to this 
system. 


- the orbital plane system. It will be assumed that the orbital 
plane precesses uniformly. 


^1 ^2 rotation Q about k^ = k 2 


^2 ^ ^3 * rotation i about i^ = i^ * 



Q is the argument of the ascending node , i is the orbital inclination. 


- the local geocentric system or desired attitude system. 

i^ i, rotation -v about jo = j/ • 

3 4 *"3 4 


Now , 


j^ is in the negative orbital normal direction 

k^ is a direction from the spacecraft to the earth*s center 
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is such that the orbital velocity has a positive component 
along i^ 

V is the true anomaly of the spacecraft. 

In practice, the desired spacecraft attitude may be that of a local 
geodetic (vertical) system rather than of a local geocentric system. For 
this case another system, whose orientation with respect to S^ depends 
upon the latitude and altitude of the spacecraft, must be introduced. 

This refinement will not be used here. 

S^ - spacecraft fixed system. This system is defined to be 
parallel to the principal axes of the main frame of the spacecraft, 

(the portion of the spacecraft not including the solar panels). The 
orientation of with respect to is determined by roll, pitch, and 
yaw as follows: 

rotation (yaw) about 

Ai A 

rotation 0 (pitch) about 

/s /s 

rotation q? (roll) about i^ = i^ . 

The nominal and desired values of , 0, and cp are zero. However, because 
of various errors in the attitude control system, these angles may vary 
between -1^ and 1^. 

S^i - the i^^ head system (see Figure 2) 

k-, -> ko . rotation B. about = jo 

7 8i -^7 -’8 

io “* rotation a. about i^ . = in . • 

1 8i 9i 

The direction will be the direction of the optical kxis of the i^^ 
head, and the angle will be called the cant angle of this head. 
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■Mf- 


PATH OF STAR 
WITH RESPECT 
TO SENSOR 

A SENSOR 

Pfield of 

\ / VIEW 



SPACECRAFT 
NADIR PATH 


Orientation of Optical Axis and Slit With Respect to S 
Two Sensor Recommended Configuration has = 02 = 62° 
= 113°, 3 , = 247°. 
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Sj^Q^j - the 1*"'^ head, slit system, 


Si ^lOij ’^°tation y. . about = j^q. 


The slit-plane of the i*"^ head has its normal. 


,-th 


The various coordinate transformations may now be written 



where 


T(i, n) 


-sin Q cos i 
-sin 0 sin i 
-cos Q 


cos fi cos i 
cos Q sin i 
-sin Q 




B(-v) 




where 


B(-v) 
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where for small cp, 0, and i|( 
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EQUATIONS OF MOTION 

A development of the equations of motion will now be given. This 
development uses, as a starting point, an expression for the angular 
momentum of the spacecraft. The writer believes that this method is 
easier to follow than that outlined by Pringle (Reference [3]) which uses 
the kinetic energy as a starting point. Assume the following: 


1. The spacecraft is composed of five rigid sections as follows: 

a. main body, 

b. three symmetric flywheels mounted such that the spin 
axis of each fl 3 rwheel has a fixed direction with respect 
to the main body, and 

c. solar array. 

2. The only factor which causes a change in the moments of inertia 
of the spacecraft is the rotation of the solar array. 

3. The center of mass of the solar array is fixed with respect to 
the main body. 


Thus, from the results of Appendix A, we can write 

4 

1 . HI. Sj) ^ H^. S^.) (2) 

i=l 


where 


H 


angular momentum of the spacecraft. 


$(T, S^) inertia tensor of the spacecraft computed with respect 
to Sy, Sy is a system with origin at the center^of mass 
of the total spacecraft. The unit vectors i^ , j^, and 
are parallel to the principal axes of the^main 
body. 

coj^y angular velocity of with respect to (the celestial 

system) . 


$(m., S .) inertia tensor of the set (m.) computed with respect to 
i mi n 1 . . . i . « 1 

Snii, S . being a system with a origin at the center of 

mass of m^. 
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i = 1, 2, 3 denotes the three flywheels, 
denotes the solar array. 

^7mi velocity of with respect to S^. 


Let 


$ (m^, 


^ ^"7 mi ^-1 ’ 

ml /mi 1 mi 


i = 1, 2, 3 


where 

d. magnitude, with respect to S- , of the angular momentum of 
^ the i^^ flywheel. 

i . spin direction of the flywheel, 
mi 

• t"Vi 

Now, let the direction of the spin axis of the i*'^ flywheel be resolved 
along the principal axes the main body. So, 


i , = b- . i_ + b^ . i- + b. . , 

mi li 7 2i *^7 3i 7 


2 2 2 

where b.. are given and b.- + b.« + b.^ 
Ji Jl J2 j3 

so that b, . 1, b..j:^0, i^i. Let 

11 ’ ij ^ 


1. The flywheels are mounted 


3 

D. = y d. b . . . 
J ^ 1 Jl 

i=l 


Hence , 


3 

I 

i=l 


S .) 
mi 


0) 


7mi 



(3) 
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Here are the components of the angular momentum of the three flywheels 
in the system. These components are known (with errors) since the 
direction, moment of inertia, and spin rate of the flywheels are known. 


The inertia tensor of the total system may be written 



where 

computed with respect to axes parallel to S^, 
mass of the solar array, 

H distance between solar array center of mass and origin of , and 

Aq , Bq , Q principal moments of inertia of the total spacecraft 
without the solar array. 


The first term in Equation (4) arises by consideration of the spacecraft 
body without the solar array. The moments A^ , , and are constant 

by Assumption 2. The second term follows from the parallel axes theorem 
(Figure IB) . Finally, the third term is simply the inertia tensor of 
the solar array computed with the origin at its center of mass. 


Hence, Equation (2) becomes 


H-HT, S,) Wjj + 


here H, KT.S^) . Dj 



2 ’ 


and mus t be resolved in 


S 


r 


’figure 1 of Appendix B. 
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EQUATIONS OF MOTION 


In order to equate the time derivative of H and the external 
torques, it is convenient to resolve H in the inertial system S^. This 
transformation has been given as 



So, A*H yields the spacecraft angular momentum components in S^^. Hence, 


T = ^ (A'H) = A'H + A'H , 
i at 

where T^^ is the sum of the external torques resolved in S^. 
T^ = AT^ = H + AA*H , 

where T^ denotes the external torques resolved in S^. Now, 
H = KT, sp (-17 -7m4^ + ^7^ (-17 -7m4^ 


( 6 ) 
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EQUATIONS OF MOTION 


where 


^17 ^7 ^2 "^7 ^3 ^7 


Given f(T, S^) , T^, I>^, and D^, the Equations (6), ( 7 ), and (8) 

determine a sixth order system of differential equations in the 
variables U)2 5 • 


In Appendix B, it is shown that 

/ 



f ^11 

- 9 i 2 



$ = 
s 

- 9 i 2 

922 


-923 


Ul 3 

■'f’23 


933 

with 





9 ii = 

Cj^ + q 

^ cos 

2p 

> 

cPi2 = 

q^ cos 

B , 



^13 = 

q^^ sin 

2B, 



CP22 

^2 ’ 




CP23 

-^2 P ’ 



CP33 

Cl - ^ 

1 

2B 

9 

where c^, ^ 

^2 ^ ^ 1 ^ 

and q^ 

are 

placement i 

of the 

solar 

array 


Thus , 


$(T, S^) 


/ cos 2p -q^ cos P 

I -q^ cos p 

\ -q^ sin 2p q^ sin p 


-q^^ sin 2p 
q2 sin 3 
- qj^ cos 



( 9 ) 
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where 


Ao + c^, 

Bo + , 

+ c^, and 


$(T, 





sin 2g 
sin p 
cos 2p 


^2 P 
0 

q 2 cos p 


2q^ cos 
q 2 cos 
2q^ sin 



( 10 ) 


In the derivation of Equation (9) (Appendix B) , it is assumed that the 
axis of rotation of the solar array is parallel to (pitch axis). 

Hence , 


‘“TmA ■ ^ 



( 11 ) 


Substituting Equations (7), (8), (9), (10), and (11) into Equation (6) , 
one obtains 
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EQUATIONS OF MOTION 


where 


= 



0 D3 -D3 


0 

-D, 



^2 = 


- 2 q, sin 2 p 


2q2 sin p ” ^2' 


0 0 
• 2 qj^ cos 2 p + 2q2 cos p 


2 q^ sin 2 p 


Hence , 



= $'^(T, S^) 


T^ - (V^ + PV 2 ) («2 



/H2 sin P 


q2 cos p 


-q2 cos p\ u^(-qj^u)2 sin 2^+q^(U2 cos p) + (Cj^-qj^ cos 2p-Bp(U2(«3‘''H2 sin P(a>2"(«3) 


q^ sin p 



2 2 , 


2 2 , 


, (-q2«^3 p " ^ 2^1 P)+(Aj^-Cj^+ 2 qj^ cos 2p)u)j^u)2-q 2 p(u)^-u)P 


2 2 , 


3 ^'^2*^1 ^ ^ 2 p)+(Bj^-Aj^-qj^ cos 2p)u)j^cu2- ^3 P (uJj^-a)2)/J 


( 12 ) 


It now remains to find the differential equations relating yaw, pitch, 
and roll to This relationship may be found by using Equation ( 8 ) or 

an addition theorem for angular velocities. Using the second method, one 
can write 


U)i7 = n - V + If ks + 0 + cp i? 
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EQUATIONS OF MOTION 



Q cos V sin i 
C(cp, 0, I -Q COS i - V 

-n sin V sin i 


Hence , 


^ sin 9 + cp \ 

^ sin cp cos 0 + 0 cos cp 
^ cos CD cos 0-0 sin cp 


V 


j 


\jf = sec 6 (o)2 sin cp + u)^ cos cp) + v tan 0 sin \)f 
+ Q(-0 cos V + sin v) sin i, 

0 = 0^2 cos cp - 0 )^ sin cp + v cos 

+ n(^ cos V sin i + cos i) , 


cp = + tan 0 (o )2 sin c() + cos cp) + v sec 0 sin ^ 

+ Q(\j/ cos i - cos V sin i) . 


( 13 ) 
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EQUATIONS OF MOTION 


Since Q is small l^/day), a small angle approximation is used in 
terms containing Q as a factor. 

The equations (12) and (13) are a system of sixth order differential 
whose solution yie Ids yaw, pitch, and roll as a function of time; given 
the constants c^, c^ j and q^, and the functions of time T^ , p, 

D^, D^, and D^. 

The flywheel angular momenta (which determine D^, D^) are not 

given functions of time, but are functions of the measured attitude and 
the control law. As stated in the introduction, an extensive effort has 
not been undertaken here to reproduce the actual mechanized attitude 
control. The control law used here will now be stated. 

It is assumed that somehow roll, pitch, and yaw and their rates are 
measured. These measured values are used to determine the momenta of 

the flywheels so that the spacecraft's attitude is nearly that of zero 
roll, pitch, and yaw. 

Since the angular velocity of (the solar array system) with 

respect to (the spacecraft system) can be written 

“7m4 " “^lm4 ’ ^17 ’ 

and since the solar array is driven so that it is nearly fixed with 
respect to the inertia system, S^; 

^7m4 ~^n ' 

So, Equation (5) becomes 
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EQUATIONS OF MOTION 


Hence, this approximation reduces Equation (12) to 

. . 2 

*'■ ^ ’ 

• 2 

Bo(i)2 = + D^tuj^ - Dj^u)3 - D 2 + (C;, + MgH - Ao)u)j (03 , 

(C„ + = T3 - D2 u,^ + + (A, - ^)u)^(U2 


where 


^7 ■ '^1 ‘7 * ^2 J7 ^3 "7 


Now, if the spacecraft is nearly correctly controlled, 
002^“^ (orbital rate) . Moreover, ^ « 6 cp 0. Hence, the 
approximate equations of motion become 

• * 2 * 

Aptui = Ti + D3V - Di - (^, - Co - MgH )vu)3 , 

BoJ >2 = - ^2 ^2 ’ 

(C;, + = T3 - D^v - D3 - (Ao - Ei^)vuj^ . 

\|f = -V Cp + (D^ 5 
0 = U)2 + V , 

cp = 0)^ + v\j; . 


Now, let 

= -(B„ - C, + M^H^) + A, . 


D3 = -(Ao - - (Co + M^H^)vcp^ + T3 , 
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EQUATIONS OF MOTION 


where cp , (u„ , and u), are determined from the measured attitude. Now, 
'm Jm 


'Im 


= (Bo - t; - )0(o)3^ - (i)J - A.; «lf^ + V D3 - + T^. 


So, 


Ao9 


-Ti + 


vD, + T, 


Let 


■^1 


= A (a cpm 



Hence , 

*’ * 1 • 

cp + bcp + acp ?=«7(T-,+ vD ) 

^ Aq i 3 

or 


cp + b cp + a 9 


2 

Ao 


(Tl 


+ ;D3> . 


So, if 4a > b , a damped motion of roll can be expected. The term vD^ 
can be nearly removed from the right-hand side by adding a "feed back" 
term which depends upon a measured flywheel momentum. This was not 
done here. In summary, the equations of motion are as follows: 


(1) Given the parameters of the solar array (Appendix B) and the 
principal moments of inertia of the spacecraft without the array 
determine the inertia tensor of the spacecraft as follows: 


= f-i n -i ' 
'1 L 3 ■ 2 


.2 X ^ .2 . 1 ,2 . 

in o)+l 24 “^2 ^ 


cos 


“]“s’ 


1.2 ,2 . 1 , 2 .. ^ 

^2 ~ ^ ^ 4 ^s 
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1 / .2 .2 ^ 1 . 2 . „ 


q2 = sin a (-j cos a + -| l) , 


Aj^ = Ao + , 


B = B, + c 


2 ’ 


:^ = Q, +M^H + c^. 


$(T, Sp = 


+ qj^ cos 20 

-q2 cos 0 

-qj^ sin 20 

-q2 cos 0 


q2 sin 0 

-qj^ sin 20 

q2 sin 0 

Cl - cos 


( 2 ) 


V = 


D3 -D2 


-D 3 0 

D 2 -D^ 



V. = 


2q2 sin p “ ^2' 

0 0 

-2qj^ cos 20 + C2 2q2 cos 0 2q^ sin 20 



E = 


2 2 , 


)l(-qiUU2 sin 20 + q^uj^ cos 0) + (Cj^ - q^^ cos 20 - 0^0)20:3 + q2 sin 0(u)2-uj3) 

2 2 

D2(-q2‘^3 sin 0)+(A3-Cj^ + 2q^ cos 20)o) j^0O3-qj^ sin 20(^3 - uj^) 

2 2 

J3(q2(JUi sin 0 + q^uj^ sin 20) + (Bj^-Aj^-qj^ cos 20)u)j^u)2-q2 ^ 
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EQUATIONS OF MOTION 


(3) The ninth order system of differential equations can now be written as 



\(f = sin © (o )2 sin cp + cos cp) + v tan 0 sin \jf 
+ Q(9 cos V - sin v) sin i, 

9 = u )2 9 " ^3 sin cp + V cos 

- Q(^ cos V sin i + cos i) , 

cp == u)^ + tan Q(uJ 2 cp + cos cp) + v sec 0 sin \jr 
+ Q(- ^ sin i + cos v sin i) , 

= 4 (a cp^ + b + (A, - + C; + , 

D„ = Bg (a 0 + b 0 ) , 

2 ° m m 

= (c; + M^H^) (a + b i J - (A, - B, + q, + M^H^) v . 

where a and b are constants; the initial conditions of (Uj^, (b^ » u)^* t 

cp, D^, are given; the torque, is a given function of i|f , 0, cp, and 

time (Appendix D) . 
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SENSOR LOCATION ABOARD THE SPACECRAFT 

An initial concern is the effect of extraneous background light upon 
the detectability of the stars whose magnitudes are near the design 
limit. The extraneous light has four sources: 

(1) direct sunlight, 

(2) earth reflected sunlight, 

(3) spacecraft reflected sunlight, and 

(4) moon reflected sunlight. 

The planets do not contribute to the extraneous background, for indeed, 
they may be used as targets. In fact, where we speak of “stars**, the 
brighter planets are included. 

The influence of the first three sources of extraneous light may 
be minimized by a judicious choice of the sensor location and viewing 
direction. This minimization can be accomplished by using the following 
constraints : 

(1) roll, pitch, and yaw are maintained at nearly zero, 

(2) the orbit is sun- synchronous , and 

(3) the solar array is the only portion of the spacecraft which 
might reflect sunlight into the sensor. 

Figure 3 is a schematic of the spacecraft which will be used to 
define the geometry of the problem. It is assumed that the roll, pitch, 
and yaw are zero so that coincides with the desired attitude S^. It 
is also assumed that the pitch axis, j^, is the axis about which the 
solar array rotates. 

Let the sensor have its origin at the terminus of Q and let Q be 
defined with respect to by 

Q = a cos Y iy ^ sin y + D k^. 


30 



SENSOR LOCATION ABOARD THE SPACECRAFT 



0 < /3 < 360“ , TO AVOID EARTH -SHINE -20" < € < -90" 
TO AVOID SUN 0 < X < 180" , D » 6* , -£ « 2 '» a , -C ,« 3* 

Ig® e\ cr = 30" 

Figure 3: Geometry of Solar Panels, Optic Origin, and Optic 

Axis Direction. 
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SENSOR LOCATION ABOARD THE SPACECRAFT 


where 

a radius of sensor mounting ring (sa 2*), 

Y azimuth of the mounting position, measured about the yaw axis, 

D distance of mounting ring (which is perpendicular to the yaw 
axis) from solar array axis of rotation, 

Because the orbit is "9:30 a.m." sun- synchronous on the average, 

n = 152.5^ + right ascension of sun. 

The direction to the sun makes an angle of approximately 60° with the 
orbital normal. The sun*s direction completes one rotation with 
respect to S^ in an orbital period. Hence, to avoid direct sunlight, 
the sensor must be positioned so that 0 < y < 180°. 

Now, let 0 denote the direction of the optical axis of the sensor, 
and let 


^ ^ 

0 = cos 6 cos e iy + sin 6 cos e jy + sin e ky , 

where 

e elevation of optical axis with respect to iy, jy plane 

6 azimuth, measured about the yaw axis, of the optical axis. 

As e increases, the sensor will point more toward the earth; and 
as e decreases, the sensor will point more toward the panels. The 
elevation of the horizon is 25° 30*. Now, to avoid extraneous light 
reflected from the sunlit earth, the angle from the optical axis to 
any ray which terminates on the sunlight earth must be greater than 
45°. This restriction is made so that a light shield of ''reasonable 
size" can be designed for the sensor. The optical axis elevation is 
set at e = -20°, 
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SENSOR LOCATION ABOARD THE SPACECRAFT 


The problem which arises because of light reflected from the panels 
must now be considered. Let P be a vector ffom the origin of S^ (on 
the axis of panel rotation) to the closest corner of the panel (only 
one panel need be considered) . Then, 


P = 


^2 

— sin p 


iy + cos 3 


A 



where 


P panel rotation angle, 0 < g < 360^, 
i, distance along to panel « 2 ' , 
panel length « 8'. 

Now let b be a unit vector in the direction of the panel's edge. 

Then 


b = sin a cos g i.^ + cos a + sin a sin g , 

where 

a panel cant angle ^ 30^. 


A vector from the origin of the optical system to a generic point on the 
panel *s edge is then 


u+pb, 0<p<>t2 » 


where 


u = Q . p = (— sin g + a cos y)iy sin y “ ^ ~ cos g) 

The cosine of the angle between the optical axis and a ray to a 
generic point on the panel's edge is 'thus, 


0-u+p0*b 

0 • d = ^ A ■ 

u + p b 


f (p) • 
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SENSOR LOCATION ABOARD THE SPACE' 


To maximize f(p) (that is, minimum angle), choose 
p = p if0<p<^2» 

= 0 if p < 0 , 

^2 P ^ “^2 * 

where 



(0 . b) (u • b) - u • 6 


0 • b = cos 6 cos p cos e sin a + sin 6 cos e cos a + sin e sin a sin p, 

2 2 2 1 2 2 

u = >t + D Y” 2 a^sinY"I^>t 2 ^ ’ 

— ^ 1 

u • 0 = — ^ 3 ^ 3 e) + a cos e cos(e - y) 

+Dsinc-^sin6 cos e > 

u * b = a (cos Y sin a cos p + sin y cos ct) + D sin a sin p - Ji cos a* 

It is difficult to proceed farther without the aid' of numerical 
techniques. We now set 

e =-20°, 

D = 6* , 

1 = a = 2 * , 

= 3 ’ , 

^2 = 

a =30°, 

and consider all values of p, 5, and y* 

/V y\ 

For any fixed value of 6 and y> the 3 which maximizes 0 • d may be 
found. This calculation yields the angle of closest approach to the panel 
envelope, that is, the surface formed by rotating the panel. 
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Figure 4 is a plot of the minimum angle between the optical axis 
and the solar panel envelope as a function of y 6. Note that to 
maximize this minimum, small y and 6 must be chosen. However, if 6 
is chosen too small, direct sunlight becomes a problem. Hence, 5 = 30° 
is recommended. With this choice of 6, it may be noted from Figure 4 
that the minimum angle to the envelope is a weak function of y. If 
£ Y ^ 20°, then the minimum angle is between 38° and 40°. Suppose 
we choose y = 0 (sensor origin in the roll-yaw plane) , then the 
minimum angle is 38°. 

Figure 5 is a plot of minimum angle between the optical axis and 
panel as a function of g. In this figure, two sensors are assumed, 
each symmetrically located with respect to the yaw-pitch plane. 

Figure 6 is a representation of the viewing geometry as seen from 
(zero roll, pitch, and yaw are assumed). The location of the optical 
axis and field of view for each sensor is shown. 

In summary, it is possible to locate two sensors on the spacecraft 
so that the angle from the optical axis of either sensor to the sun 
is greater than 50°, the angle to the horizon is greater than 46°, and 
the angle to the panel is greater than 38°. It is felt that an 
extraneous light shield of not unduly large size can be designed with 
these constraints. The angles and which define the position of 
the i^^ optical axis (these angles were defined in the section on coordinate 
systems) are 


= 113 ° 
P2 = 247° 



azimuth of i optical axis measured with respect to the pitch- 
yaw plane 

a. coelevation of optical axis with respect to the roll-yaw 
plane (Figure 2) 
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minimum equator maximum spacecraft minimum 

LATITUDE I/=I80» LATITUDE AT EQUATOR LATITUDE 

1^ = 90® ^t/ = 270® 4i/= 0® l/= 90® 


SENSOR LOCATION ABOARD THE SPACECRAFT 


o 
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Figure 5: Minimum Angle Between Optical Axes of Each Head and Solar Panel as a Function of 

Solar Panel Angle. 




HEAD* 2 
FIELD OF VIEW 


SENSOR LOCATION ABOARD THE SPACECRAFT 
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Figure 6: 


Viewing Geometry as seen from S^. 
-Icy is the Zenith Direction. 


i^ Roll, Pitch, Yaw. 
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EQUATIONS IMPLIED BY A MEASURED TRANSIT TIME 


Recall that the system S,„.. is a coordinate system associated 
th th 

with the j slit of the i sensor, and the slit is in the ^lOij^ 

plane. Hence, the instant a star is in the slit plane (that is, at a 
transit tim^, the direction to the star is given by 


Uio = cos T1 jg. + sin T1 . . 

The direction to the same star may be written as 


= cos 6 cos a i^ + cos 6 sin a + sin 6 j 


where o' and 6 are the right ascension and declination, respectively, 
of the transited star. Hence, from Equation (1) , 



B(Vij)S(ai,Bi)C(cp,0,\^)B(-v)T(i,n) 



(15) 


The parameter, T], may be eliminated from the above systems by simply 
considering the first of the equations. Thus, 


0 = B^(v.j)S(a^, 3 pC(cp, 04 )B(-v)h , 


(16) 


where 


= first row of B , 


: 

h =1 j = T(i,Q) 

\ ^3 / 
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EQUATIONS IMPLIED BY A MEASURED TRANSIT TIME 


That is, h^, h2 > and are the components of the star in S^. 

Equation (16) represents one equation in the three unknowns cp(t), 
0(t), and \l/(t), where t is the transit time. This equation is obtained 
from each transit of a known star. If no further information is given, 
each new equation simply adds one equation and three unknown angles - 
the roll, pitch, and yaw evaluated at each isolated transit time. 

The physics which governs the motion of the spacecraft must be 
invoked to develop a time-dependent characterization of the attitude 
which involves just a few unknown parameters. This ’’attitude model” 
may then be used to internally couple roll, pitch, and yaw. This 
important point will now be considered in some detail. 
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THE ATTITUDE MODEL 

Since attitude of the spacecraft as a function of time is completely 
determined by the spacecraft inertia tensor, flywheel angular momenta, 
external torque, and initial attitude, it is natural to attempt to use 
a set of simplified equations of motion to supply a model. Simplified 
equations rather than the more exact equations are necessary since the 
exact equation would require excessive computer running time to reduce 
operational data. 


One simplification which seems necessary is that the solar array 
remains fixed with respect to an inertial system. This approximation 
yields the system of Equations (14). A comparison of the solution implied 
by this approximation and the true solution is given in Figure 8. 

Thus, the approximation produces an error of .03 degrees after 
50 seconds. Such an error is on the border line of being accep- 
table . 


To use the approximate equations as a model, additional assumptions 
must be made. The external torque, flywheel angular momenta, and flywheel 
torques must be somehow supplied. The flywheel angular momenta and 
torques can be measured (with errors) . The external torque may be 
obtained^ by assuming it constant in the data gathering interval. Such 


The external torque may also be computed by assuming Equations (14) 
with roll, pitch, and yaw, on the average, zero. This yields 

(t2 - tp <T^> = <D^> - ; <D3> , 

(t2 - t^) <T2> = <D2> . 

<"^2 - <^3> = <°3> + ^ <°1> ’ 


where 


(t2 


<Ti> average value of T^ in [t^, , 

tp <i).> = D.(t2> - , 

<°i> average value of in 


The difficulty with this method of estimating the torque is that it tends 
to produce poor estimates if the external torque were to change suddenly 
or often. 
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IMPLIED BY APPOXIMATE 
EQUATIONS AND MEASURED FLYWHEEL 


THE ATTITUDE MODEL 



Figure 8: Yaw, and Two Approximations of Yaw as a Function of Time. 




THE ATTITUDE MODEL 


a model then contains 9 unknowns; so a minimum of 9 transit times must be 
held in the data reduction filter. To obtain a reliable attitude in the 
presence of typical transit time errors, approximately 20 transits must 
be held in the numerical filter. Later (pages 57-66), it will be shown 
that for the proposed instruments; transits are gathered at a rate of 
approximately 0.3 transit/sec. Thus, on the average, a data gathering 
interval of 60 seconds must be used to obtain 20 transits. 

Figure 8 is a plot of yaw as implied by a solution of Equations (14) , 
in which the flywheel momenta and torque are measured (107o RMS error) and 
the external torque is set equal to its average value. 

Since a model which is based upon a set of simplified equations of 
motion implied a somewhat complicated data reduction scheme and may 
produce border line accuracies, such a model was abondoned. In retrospect, 
perhaps this was premature. 

Constant Yaw, Pitch, and Roll Model 

Instead of assuming a model which allows for a complicated attitude 
motion, an investigation was made as to the results attainable from a 
model which assumes yaw, pitch, and roll are constant during the data 
gathering interval. It turns out that if certain refinements are made, 
the constant model does generally provide a background from which quite 
accurate attitude estimates can be made. This general point will now be 
discussed. 

Good Geometry 

Assume an attitude estimate is required at some given time, T^. Since 

only three transits are necessary to determine yaw, pitch, and roll, it 

seems natural to use the three measured transit times nearest to T . These 

P 

three transit times will generally span an interval of about 10 seconds; 
so if the attitude is not changing rapidly, the constant model cannot 
introduce a large error. 

A major difficulty with this approach, however, is that the geometry 
implied by the three measurements may be ill-conditioned. The nature 
and development of this ill-conditioning will now be discussed. 
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THE ATTITUDE MODEL 


Consider the instant a star is in a slit. Three orthonorraal 

A 

directions may now be defined (Figure 9): s, the direction to the 

/V ^ ^ yv A 

Star; r, the slit normal; and p = r x s. The direction p will be 
called the "preferred direction". Note that the spacecraft may be 
rotated with an arbitrarily large angle about s or r without causing 
the star to leave the slit. However, a rotation about p causes s to 
move in a direction perpendicular to the slit. Each transit thus 
establishes a preferred direction about which attitude motion is not 
allowed at the transit time. If only three relatively close transit 
times are measured, the best geometry would be such that the three 
preferred axis form an orthonormal triad. If the three preferred 
axis lie in a common plane, then the geometry is ill-contained. 

Suppose a single head is used, then each preferred axis is 
approximately 90^ from the optical axis (Figure 10). Hence, no 
preferred axis has a large component about the optical axis. This 
in turn implies that the single-head sensor yields a geometry which 
is inherently poorly conditioned since attitude about the optical 
axis is poorly defined. 

If only three stars are used and the half-f ield-of-view angle is 
f, the placement of the stars which yields the optimum geometry is 
given in Figure 11. The optimum geometry is that which yields a maximum 
volume parallelepiped whose coterminal edges are the preferred axes. 

For f = 10° (the recommended ha If- fie Id) the optimum configuration 
yields a parallelpiped of volume 0 . 202 . 

Now consider a two-headed sensor, for simplicity, and assume the 
angle between the optical axes is 90° (the recommended configuration. 

Figure 6 , has a separation of 108° 44*)- The configuration S^, > and 

So in Figure 12 will yield a well-conditioned geometry since the respective 

/V /% yv /s. 

preferred axes of this triad are approximately ^ ^ 2 ’ ^ 2 ’ *^ 1 ' 
configuration S^, and S^, however, will yield an ill-conditioned 

A A 

geometry since its respective preferred axes are approximately 0 .. x 0 ^ , 

O 2 , 0^ X O 2 . Rotation about 0^ is thus poorly determined by the 
configuration 82 ? and S^. 
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A 

S STAR DIRECTION 



A 

P, PREFERRED DIRECTION 

Figure 9: The Preferred Direction, 


OPTICAL AXIS 



Figure 10: The Preferred Axes for the Single-Head Sensor. Note 

that Each Preferred Axis Has a Small Component on the 
Optic Axis; Hence, Attitude About the Optical Axis is 
Poorly Determined. 
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THE ATTITUDE MODEL 


For the recommended two-head sensor the optimum configuration of 
three stars yields a parallelepiped of volume 1.0 (compared to .202 for 
the one-head sensor). Hence, the two-head sensor has a significantly 
higher probability of providing a well-conditioned geometry with a small 
number of stars than does the single-head sensor. 


In practice, all transits are not gathered simultaneously as was 
assumed in this discussion. Instead, they are typically gathered over 
a minute. This effect increases the directional region from which transits 
can be obtained and thus increases the probability of obtaining a good 
geometry. 


The general method of choosing, from a given set {^j] of ordered times, 
the transit times used to compute the attitude at a given time, T , is as 
follows : 

(1) Choose three transit times, t., t . . - , and t, such that 

1 L+1 k 


t , = [t I min (T - t .) ], 
1 ' ' t.<T P J • 


J P 


t = {t 1 min |T - t \]. 

^ i+1 P J 


(2) If these three times yield a good geometry, then the attitude 
is computed with these times. If they do not yield a good 
geometry, the next closest time to T is chosen. For example. 


t = 



min 

j?^i, i+1, k 



The geometry with these four times is investigated, and the 
process continued. 


A criterion as to whether or not a given set of transit times produces 
a good geometry can be obtained by examining the coefficient matrix of the 
linear system of equations, or by examining the star-slit position at each 
transit . 


Weighting of Measurements 

If the constant yaw, pitch, and roll models were accepted and the 
variance of each measurement transit time were known, the Guass-Markof f 
theory from statistics (Reference [4]) gives a "weighting" of each equation. 
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It is known, however, that yaw, pitch, and roll are not constant; hence, 
those measurements which are nearer to the time at which the attitude is 
sought must be trusted more than those far removed. Unfortunately, 
there is no simple theory which gives a method of selecting the weights 
for the present problem. Instead, the weights were selected empirically 

Let T be a time at which an attitude estimate is desired. Also, 

P 

let T. = t. - T , where t. is measured transit time. Let 

ip 1 p L 






- 1 


10 


ip> 


A graph of f(T) is given in Figure 13. It was felt that f(0 was 

undesirable to use as a correction to the weights since the first 

moment about T may not be zero. A correction was therefore added so 
P 

that the first moment is zero, and the sum of the squares of the 
correction is a minimum. Thus, the weights are determined by 


w . 
iP 



T . . 


Li ip 

i=l 




P 2 

ip 


2 

L 


i=l 




where n^ is the number of measurements used to yield an estimate at 
Tp. It turned out that using the first-moment correction, on the 
average, did supply a slightly improved estimate. 

Estimation of From Flywheel Outputs 

For the nominal case, each flywheel rotates at an almost constant 
rate. The slow variations in their rates are caused by slow variations 
in the external torque. However, it may happen that the horizon 
sensor may supply an erroneous measurement due to a variation in the 


49 



o w 


THE ATTITUDE MODEL 





CO 

0) 

u 

cd 

(U 


c 

<u 

u 

o 

CO 


<u 

CO 

o 

H T3 
Q) 

0) U 
N 'rH 
•I-I CO 
1-H <y 
•H Q 
4-i 

P CO 
<U 

M (U 
O TP 

S o 

U 
O -H 
U U 
AJ 

CIO < 

c 

•H <U 

4J 

W) 

43 
(U O 

4J 

TP Cd 

o 

rC Q) 

4J 6 
Q) 'H 

S H 


CO 


(U 

}-l 

3 

bO 


50 



THE ATTITUDE MODEL 


earth’s radiation output. These erroneous measurements will cause a more 
rapid change in the flywheel rotation rates, and hence in the angular 
velocity of the spacecraft fixed system, S^, with respect to the inertia 
system, S^. A method of estimating this angular velocity, from 

the time history of the flywheel output will now be given. 

Assume the following: 

(1) The solar panels are always at a fixed position with respect 
to an inertia system, S^; 

( 2 ) 0 )^ ^ ^ Oj and -v (orbital rate) ; where and a)^ 

are the components of with respect to S^; 

(3) Roll, pitch, and yaw are 0. 

Then, approximate equations of attitude motion were found to be (Equation (14)) ^ 
Ao(i)i = Ti + D^v ' ^ (®o ~ Cq - M^H^)vu)2 " 

= -D 2 + T^, (17) 

(C, + = T 3 - D^v - D 3 - (A, - 

\j/ = -VCp + 0)^ , 

0 = U )2 V > (18) 

Cp = 0)^ + > 

where A^ , Bq , and are the principal moments of inertia without the 
panels, and 

mass of the panels, 

H distance between center of mass of the spacecraft with 

panels and the center of the mass of the panels, 
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^ 1 *^ 2 ’*^3 
•il, Q, cp 


angular momentum components of the flywheels, measured 
with respect to the spacecraft system, , 

components of the external torque resolved in , 
yaw, pitch, and roll, respectively. 


First of all, the external torque may be estimated by assuming the 
following: 


(1) ^ ^ where <^x> is the average 

value of X in [t^, t^]; 

(2) T^, and are slowly varing functions of time. Thus, 

<ij> (tj - tj) = <i^> -;<)3> . 

< C ^ 2 ^ (^2 " " ^^ 2 ^’ 

<'^3> <'=2 • "l> ' <®3> ■ 

where 

<4 - <“l> ■ ”l<'=2> - 


Now, suppose there exists an interval of time < t < T 2 such 
that for t in this interval 


Ti + D 3 V 


- D, 






0 . 


Then, we may assume = u)^ ” CO 2 “ this interval. Let t > X 2 > 

then the solution of the differential Equations (17) is 
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(O 2 (t) 


■ill - D^Ct)] , 


where 


$(t - s) = 



cosh X(t - s) 


— sinh X(t - s) 
a 


7* sinh X(t - s) 

D 


cosh X(t - s) 


A,f^(s) = T^(s) + vD^(s) - D^(s) 


(Q + )f3(s) = T^Cs) - vD^(s) - B^(s) 

2 


Bb - C, - M H 
3 > ^ .33 , 


b = ^ ~ « .29 . 

Q + M H 
° s 


X = V \/ab « .31v 3.1 X lO"'^ rad/ sec . 


If t - T 2 < 300 seconds, further approximations are possible so 


that 


oo^(t:) = 


f^(s)ds , 


cu 


I^(t) = J f2(s)ds. 

■^2 
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The yaw, pitch, and roll rates may then be approximated from Equations (18). 

Because the D. cannot be measured without error and various 
1 

approximations were used, the estimates of o)^ have errors. Below, a 
graph is given which compares the true pitch rate with the estimated 
rate . 
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TIME (sec) 

Figure 14: Estimated Pitch Rate and True Pitch Rate as a Function of Time* 

1640 < t < 2400 Seconds 





DETERMINATION OF ATTITUDE 


DETERMINATION OF ATTITUDE 

In this section, the contents of the previous sections will be 
summarized and a final method of attitude determination will be given. 


The basic constraint Equation (16) may be written 
0 = D(y^j, pp C(\|;, 0, cp)q , 

where D and q are unit row and column vectors, respectively, and are 
given by 


D(Yij» pp = S(a^, pp 


q = B(-v)h 


So, for small yaw, pitch, and roll 



with 


(19) 


a, = -D q, and 

I 

\ = ^^ 1^2 - ‘^ 2 ^ 1 ’ ^^ 3^1 ■ ^3 - ^ 3*^2 > 

where the q*s and d*s are the components of q and D, respectively. 


If the components of 
from the flywheel output, 


the angular velocity, 
then 


are estimated 
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cp ^ 0)^ + 



\lf ^ " vcp • 

The solution of this linear system of differential equations is 

t 

cp(t) = co(Tp) cos \ + \lf(Tp) sin X + JCu)j^(s) cos v(t-s) + ^ 3 ( 5 ) sin v(t-s)Jds, 

T 

P 


9(t) = -v(t - Tp) + u) 2 (s)ds + 9(Tp) , 

T 

P 


(t) = \li(T ) cos X - cp(T ) sin X + r [-u), (s) sin v(t-s) + ^ 3 (s)cos v(t-s)]ds, 
P P J 


where 

X. = (t - • 

• o 

Since v lO” rad/sec and |t - < 100 sec, X < .1 rad. Hence, 

to a first order approximation, 

cp(t) «cp(Tp) + (t - Tp)(w^(Tp) + vil,(Tp)) 
e(t) ^(Tp) - v) (t - Tp)+ 0(Tp) 

\lf(t) (t- Tp) (u) 3 (Tp) - vcp(Tp)). 

Substituting these results into Equation (19), one obtains 
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/♦‘V\ 

■ ^^1 ’‘4*13’ *42’ *43 ■ '’4*41' f j 

V'P'Tp)/ 


If a number of measurements are to be processed to obtain the 
attitude at T^, a system of equations can be written as follows: 


^(Tp) 

b-B I 9(Tp) |, b(r.p X 1), B(i.p X 3). 


These equations are then weighted by the criterion discussed in the 
previous section to obtain 


Wb = WB 



where W is a diagonal matrix 


(n X n ) . 
P P 
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STELLAR AVAILABILITY AND DETECTABILITY 

At this point, only the direction of the sensor *s optical axis with 
respect to Sy has been specified. A number of other parameters must, 
however, be specified in order to discuss the probability and accuracy 
of stellar detections. These as-of-yet unspecified parameters are listed 
in Table V. Also listed in Table V are numerical values of these 
parameters. The procedure employed in choosing these numerical 
values is obtained through the iterative process shown in Figure 15. 

The total process is described in References [5], [6], and [7] and 
will not be repeated here. 

In Table V, the parameters above the double line are independent. 

Those below the double line are a function of the orbit, distribution 
of position and magnitude of the stars, and the independent parameters. 

Figures 16 and 17 give the slit pattern for the passive and active systems. 
Note that the passive system has a much smaller central dead zone than 
does the active. Also, note (Table V) that the width of the slit on 
the mask of the passive system is greater, by a factor of 5, than that 
of the active. Each of these effects occurs for the following reason: 

In order to make the transit duration independent of the star position 
along the slit, the mask for each slit is wedged-shape (Figure 16) with 
vertex at the spin axis projection upon the plane of the mask (focal 
plane) . For the passive system, the spin axis is (approximately) the 
satellite pitch axis which is 62^ from the optical axis. Hence, within 
the field of view, the sides of the slit projection are nearly parallel. 

The active system, however, is such that the spin axis and the optical 
axis are (approximately) the same. Hence, the active system must possess 
a relatively large dead zone and rotational slit width to avoid the 
convergence of the slit sides near the optical axis. 

Consider the following instrument parameters as given: aperture 

diameter, area of slits on the photodetector, solid angle subtended 
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Figure 15 : Iterative method of obtaining instrumentation 

parameters. The computer program I is an outgrowth 
of that reported in [7]. The program II is that 
discussed in the previous section. 
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by the slits, optical efficiency, photocathode sensitivity, slit geometry, 
and electrical filter. These inputs, together with the directional and 
intensity distribution of the stars allow the computation of the following 
outputs as a function of the instrument spin period: background noise at 

filter output, scanning noise* at filter output, signal noise at filter 
output, signal- to-noise ratio, standard deviation of the transit time 
error, expected number of false transits per unit time due to scanning 
and background effects, number of primary photoelectrons which arrive 
at the photodetector as a function of stellar magnitude, and probability 
of detection as a function of stellar magnitude. The method of 
computation is explained quite completely in Reference [8]. 

One output from Program I (Figure 15) is the stellar detection 
probability as shown in Figure 18. This output is obtained again by 
completion of the iteration indicated in Figure 15. A second output is 
the estimated transit time error. 

The standard deviation of the measured transit time may be found 
as follows: 


a(m, P) = cT^Cm^CP)) 


( 20 ) 


where 


magnitude of star being detected, 

magnitude of star detected at the threshold level, this 
level being a function of the instrument spin period, 

standard deviations of the measured transit time if a 
star were detected at the threshold level**. 

Values of a^(m^(P)) are given in Table VI. 


iU3,(P) 

O.J, (nVj, (P) ) 


*Noise due to dim stars which are close to the slit at the time a target 
star is detected. This noise source may be eliminated by rejection of 
such contaminated target stars. 

**The threshold level is chosen so that the probability of detecting a 
star of magnitude m^(P) is 0.5. 
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PITCH AXIS 
TERMINUS ^ 



THE EXTENSION 
OF SLIT 3 



ROTATION SLIT 
WIDTH W, = 2 mTn. 


r^=i«, rg= I0» 

ANGLE BETWEEN OPTICAL 
AXIS AND PITCH AXIS , 

(T = 62 ® 


AVERAGE TIME FOR 
STAR TO PASS FROM 
A TO B IS 365 sec. 

EXPECTED No. OF 
TRANSITS FROM A 
SINGLE STAR =4.3 


Figure 16: Slit Pattern for Passive System. If Two Heads are 

Used, Slit 5 is Continuous in Second Head. 
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Figure 17: Slit Pattern for Active System. 
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TABLE VI 

m^(P) and a^(m^(P)) 


P (Seconds) 

nr^(P) 

(tOj, (P) ) (Seconds) 

5 

3.5 

.000081 

10 

4.0 

.00012 

20 

4.5 

.00035 

30 

4.75 

.00053 

6200 

5.00 

.041 
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Figure 18: Detection Probability as a Function of Stellar Magnitude and Instrumentation Spin 

Period. P ® 5, 10, 20, 30 (2" Aperture), and 6200 Seconds (1" Aperture). 





STELIAR AVAILABILITY AND DETECTABILITY 


All orbital elements except the longitude of the ascending mode, 

Q, have been specified. Hence, if this element is given, the orbit 
is completely specified and thus stellar background is known. This 
information together with the detection probabilities given in 
Figure 18 are sufficient to compute the expected number of stellar 
detections per second. This information is given in Figure 19. 

First note (Figure 19) that the average number of stellar detections 
depends more strongly on Q as the instrument spin period is decreased. 

This occurs because the shorter periods imply detections of the brighter 
stars only (also, of course, more detections per det:ected star), but the 
brighter stars tend to be distributed in clusters. That is, the dimmer 
stars tend to be more uniformly distributed that the brighter stars. 

Next, note that 0.2 transits/sec will be obtained with the passive 
system and 0.4 transits/sec will be obtained with 30 sec spin period 
active system. A first thought might be that the 30 sec active system 
yields attitude whose errors are one-half those of the passive system. 

It will be shown that this is, in general, not a good approximation. 

This point will be discussed in detail later. 

Recall that the apertures of the passive and active systems were 
fixed at 1” and 2 ”, respectively (Table V). Thus, the active system demands 
a larger instrument. Of course, it may be argued that a smaller active 
instrument (than 2”) can provide more transits/sec than the proposed 
passive system. It seems, however, that if one is to pay for the added 
instrumentation necessary for the active system, the resulting attitude 
estimate must be better than a marginal improvement over that afforded 
by the passive system. This criterion fixes the aperture of the active 
system at 2*' and the maximum spin period of the active system at 30 seconds. 

Finally, note from Figure 19 that, except for a certain set of Q, 
little or no improvement in number of transits/sec is obtained from 
spin periods less than 10 seconds. 
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Figure 20 gives the general method used in this study to determine 
the accuracy to which the stellar transits determine the spacecraft's 
yaw, pitch, and roll. Note that four distinct attitudes are used. 


(1) The attitude as measured by the onboard attitude control 
system. This attitude estimate is used to drive the 
flywheel, and in this study it is considered a somewhat 
poor estimate of the true attitude. The parameters used 
to drive the fl 3 ^heels are thus 


9 

9 

i 

9 


m 


+ 

Allf 

m 

Q 

+ 

AS 

m 



Acp 

m 


+ 

Ai; 

m 

Q 

+ 

a6 

m 

9 

+ 

Acp 


where > AU , Acp, A\[f > A9 , and Acp are chosen (as a function 
of time) input errors. In this study, the error in control 
attitude rates were chosen as A\|f = A® = Acp = 0? and the 
errors in control attitude were chosen as the functions 


At(t) = 2C. 

A9(t) = ^ .05° sin 2vt, 

Acp(t) = , 



t 


i-1 


+ 


e < t < t^ ; 
e < t < t. ; 
e < t < t . . 


In the interval t.<t<t.+e, a smoothing formula was used 
so that the error in control attitude is continuous. Numerical 
values of the C's and t*s are given in Table VII. 

These errors were chosen to simulate constant bias errors in 
the horizon sensor. These constant bias errors change rapidly 
to simulate changes in the observed radiation output of the 
earth as the spacecraft orbits the earth. The sinusoidal term 
was added to the control pitch error to simulate an uncompensated 
error due to the earth oblateness. 


69 




RESULTS 



Figure 20: Method Used to Compute True Attitude and Star 

Sensor Attitude. 
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TABLE VII 

PARAMETERS WHICH DEFINE THE CHOSEN ERROR IN THE 
CONTROL ATTITUDE (t^ = 0) 


i 

t . (sec) 
1 

C^(deg X 100) 

1 

170 

- 8 

2 

410 

- 2 

3 

1000 

2 

4 

1670 

-10 

5 

2000 


6 

2670 

4 

7 

3330 

6 

8 

4170 

- 4 

9 

5000 

2 

10 

cx> 

0 
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(2) The true attitude. This attitude is obtained as a solution 
of the differential equations with chosen initial conditions 
and chosen input parameters as previously given. A plot of 
resulting yaw and pitch is given in Figure 21, An explana- 
tion of its behavior is in order. To do so, let us concentrate 
on yaw in the time interval 450 < t < 2000 seconds. 

In the interval 450 < t < 1000 seconds, = .04^. The control ^ 
law is designed to drive = \|; + A^ to zero; hence, ^ ^ -A^ = -.04 

as is shown in Figure 22. At t = 1670 seconds, the error in control 
yaw jumps to Aij; = -.10^. However, \|; does not immediately jump 
to 0.10^. Instead, there is considerable overshoot and a value of 
= 0.80^ is reached. After about 300 seconds, the jump in the 
control yaw error is damped out so that in the interval 
1200 < t < 2000 seconds, \jf ^ -A\|/ = 0.10^. At t = 2000 seconds, 
another jump in the control yaw error occurs. The process is then 
repeated, this time with more overshoot because of the greater jump 
in control yaw error. 

(3) The model attitude. This attitude is the functional form assumed 
for yaw, pitch, and roll over the data gathering time interval. 

This model attitude has been discussed previously. 

(4) Finally, there is the attitude as estimated by the transit time ^ 
and the flywheel outputs. This attitude is denoted by \|; , 9, and cp. 
In this study, an attitude estimate will be supplied every 20 
seconds over one orbital period. An attitude estimate could be 
supplied more often with no difficulty. 


The accuracy with which the attitude can be determined by use of the 
stellar transits will now be considered. The passive system and various spin 
rates for the active system will be investigated. For both the passive 
system and the active system, the stellar field implied by Q = 0 was chosen. 
This choice was made because Q = 0 generally gives the smallest number of 
transits per orbital period over all choices of Q (Figure 19). This choice 
thus, on the average, will yield larger attitude errors than other choices 

of n- 


Passive System 

Salient error results for the determination of attitude with the passive 
system are given in Figures 23, 24, and 25. Figure 23 gives cumulative 
distribution of the roll error. To obtain this figure, 310 cases were 
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Figure 22: True Yaw, Error in Control Yaw, Yaw Flywheel Angular Momentum, and Yaw External Torque 

for 450 < t < 2000 Seconds. Control Law is Such that ik + Ail/ is Driven to Zero. 
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Figure 23: Cumulative Distribution Function of Roll Error. Passive System, One and Two Heads, 

Aperture 1". Errors in Yaw and Pitch are Slightly Smaller than the Roll Error. 
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Figure 25: Yaw and Pitch, and Two Estimates of Yaw and Pitch for 1640 < t < 2040 Seconds 





RESULTS 


examined (an estimate every 20 seconds over the orbital period) . Note that 
the passive system yields a roll error of less than 0.02^ in 87 percent of 
the cases if a two-headed sensor is used. The corresponding number for 
the s ingle-headed sensor is 58 percent. 

Figures 24 and 25 compare two estimates of yaw and pitch with the 
true yaw and pitch in two different time intervals. Note that the attitude 
obtained with the aid of the rate estimates as described previously yield 
a significant improvement over that which simply used a constant model 
attitude. Roll estimates are not shown in Figures 24 and 25, but they are 
quite similar to the errors shown. 

Active System 

Results derived from active system parameters with spin periods of 
10, 20, and 30 seconds are summarized in Figures 26, 27, and 28. Figure 26 
shows the cumulative distribution function of the absolute value of the 
roll and yaw errors for the two-headed active system. The corresponding 
function for the pitch errors is not given, but it is quite similar to 
the plotted functions. From this figure, one notes that the errors are 
somewhat independent of the spin period. Also, the errors are quite small 
in that the probability is 0.96; that the error in any attitude angle is 
less than 0.03^. 

Figure 27 is of the same format as Figure 26 except it is obtained 
from the parameters of the single -head system. Note that there is a 
significant change as the spin period is varied. For the single-head, 
the 10 second period yie Ids larger errors than that of the 30 or 20 second 
periods. This increase in error occurs even though the data rate (Figure 19 
and Table VIII) is approximately the same as that of the 30-20 second 
period. It occurs because at the faster spin rate only the brighter stars 
can be detected; hence, there may exist relatively long intervals of time 
over which only one or two stars are detected. This state implies large 
errors and is more probable with one head than with two heads. 
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jure 26: Cumulative Distribution Functi 

Periods. Two Heads, 2" Apertu: 
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Figure 27: Cumulative Distribution Functions of Roll and Yaw Errors 30, 20, and 10 Second 

Periods. Sinble Head, 2" Aperture, Active System. 






RESULTS 


A central question which was to be answered by this study is: 

If an active instrument is to be used, what is the optimum spin rate? 

An answer may now be given. 

Suppose, for the moment, that optical aperture diameter is fixed at 
two inches and the field of view is fixed at 20^, then periods much 
longer than 30 seconds should not be considered. In this case, the 
accuracy of the system approaches that of the passive system. Why 
consider an active system if a small passive system yields a comparable 
accuracy? Periods much shorter than 10 seconds should not be considered. 
In this case, time intervals can exist in which only one or two stars are 
detected if a single-head sensor is used. But, even if the spacecraft is 
fitted with a two-headed sensor, periods will exist in which the moon 
will cause one head to be inoperative. Hence, single-head operation must 
be considered. Shorter spin periods could be considered if the aperture 
diameter were greater than two inches or the field of view were greater 
than 20°, for each of these changes would imply that more stars would be 
detectable. The first change implies a larger instrument, and the second 
implies a more difficult optical and stray- light shielding problem. 

Neither change is recommended. Our recommendation for an active system 
is one which consists of two heads, each with a spin period between 20 and 
30 seconds. 

Additional results are given in Table VIII. This table lists the 
standard deviation of each estimated attitude angle, the number of stars 
which were detected in an orbital period, and the total number of transits 
gathered in an orbital period. This information is given for the two-head 
and single-head sensor and for three active systems. 
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AREAS FOR ADDITIONAL INVESTIGATIONS 

This study may be considered as a fairly complete analysis of the 
feasibility of providing a high accuracy, operational, near real-time 
attitude determination system for the ERTS spacecraft. However, several 
areas have been uncovered during the study in which system changes may 
improve the accuracy of the system. These points and additional areas 
for study will now be commented upon. 

Ins tr umen ta t ion 

It is possible to consider different slit patterns for both the 
passive and active systems which should yield smaller attitude errors 
at the expense of more difficult optical designs. Smaller errors may 
not be of importance for the two-headed active system, for the present 
system does provide quite satisfactory results. However,. it may be 
desirable to improve the accuracy performance of the passive and single- 
head active systems. 

A possible slit pattern improvement for the active system is shown 
in Figure 29. This figure shows a three-slit configuration similar to 
that successfully used by GDC on the ATS- III spacecraft (Reference [ 9 ]). 

This configuration provides a coelevation measurement (angle between a 
star's direction and the optical axis) as well as the azimuth measurement 
obtained with the single slit. The coelevation measurement is obtained 
because the time between successive transits of the same star is a function 
of the coelevation. That is, three closely occurring transits will 
characterize a star far from the optical axis. As the star is moved toward 
the optical axis, the transits will become more widely spaced. 

This three slit configuration also makes possible a simpler star 
identification problem than that implied by a single slit. This simpli- 
cation occurs because more information is obtained concerning each star's 
direction so that each candidate star is required to pass a more severe test 
before it is paired to the transited star. The disadvantage of the three- 
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PRESENT ACTIVE 
SLIT PATTERN 


Figure 29: Possible Improvement 

System. 



PRESENT PASSIVE 
SLIT PATTERN 


AREAS FOR ADDITIONAL INVESTIGATIONS 



POSSIBLE IMPROVEMENT 


of Slit Pattern for the Active 



POSSIBLE IMPROVEMENT 


Figure 30: Possible Improvement of Slit Pattern for the Passive 



AREAS FOR ADDITIONAL INVESTIGATIONS 


slit configuration is that the optical system must possess small radial 
distortion as well as small tangential distortion. For the single-slit 
configuration, radial distortion does not give rise to errors. Also, 
the background noise is greater for the three-slit configuration than 
the single slit configuration. 

Figure 30 gives the present passive slit pattern and a possible 
improvement. The suggested improvement will provide more transits per 
star and thus a more accurate system. Moreover, it could ease the star 
identification problem. This can be done by designing the pattern so 
that the time between successive transits of the same star is constant if 
yaw, pitch, and roll are zero. Again, the optical design problems are 
greater with the suggested change. 

Star Identification 

The star identification problem has not been considered in this study. 
Each slit pattern presents a separate problem. It is most difficult for the 
present passive slit pattern, for in this case the star as well as the slit 
must be identified. The problem may be considered solved for the suggested 
change to the active slit pattern. 

Analysis 

It may be possible to improve system performance in the presence of 
low data and rapid changes in attitude. This could be done by use of 
approximate equations of motion and a Kalman- type filter instead of the 
model and filter used in the present study. The disadvantage of this 
suggested approach is that systems of differential equations must be 
solved numerically in the data reduction process. Thus, a longer computer 
running time will be required. Perhaps this more time-consuming reduction 
should be used only if a low data rate is coupled with a rapid change in 
attitude. 
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CONCLUDING REMARKS 


Ins trumenta t ion 

Two classes of instrumentation were considered, a passive system 
and an active system. The passive system consists of five radial 
slits which are fixed with respect to the spacecraft. Stars transit 
the slits due to the slow orbital motion of the spacecraft. For this 
system, a one-inch diameter optical aperture can be used and fifth 
magnitude stars detected. The system is very attractive from the 
instrumentation view point since a small, light weight, low power 
consumption, reliable instrumentation with no moving parts is possible. 

The active system consists of a single radial slit (three times 
wider than each slit of the passive system) which is made to rotate 
about the sensor optical axis. This rotation is produced by simply 
rotating the small thin reticle upon which the slit is photoetched. 

Thus, the active system gathers stellar transits at a greater average 
rate than the passive system. However, its optical aperture must be 
larger (two- inch diameter) and a method must be provided to measure the 
orientation of the slit about its rotation axis (angle encoder or 
constant drive motor) . The recommended spin period of the reticle is 
20 to 30 seconds. At these periods, stars of magnitude 4.5 to 4.75 will 
be reliably detected. 

For both systems, the recommended detector is a photomultiplier. 
This recommendation is made simply to make possible a smaller optical 
aperture than if a solid state detector were used. 

Two sensor heads are recommended for both the passive and active 
systems. Two heads provide a stronger measurement geometry, greater 
data rate, redundancy, and system operation in the presence of the moon 
(the moon can saturate only one of the two heads at a time). The field 
of view of each head is 20^ with a central dead zone of 3^ for the 
active system and 1^ for the passive system. 
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To yield maximum angular distance from the earth (which may be 
sunlit) , sun, and solar array, the following head placement should 
be used. Each head should be placed 62^ from the pitch axis. The 
projection of the first head's optical axis on the yaw-roll plane is 
113^; the second head's, 247^ (see Figure 2). For this case, the 
angle from the optical axis of either sensor to the sun is greater 
than 50°, the angle to the horizon is greater than 46°, and the angle 
to the panel is greater than 38°. 

Data Reduction 

A batch processing, least squares estimator is recommended. The 
model attitude is of the form 

cp(t) = r^(Tp) 4 (t - Tp) + vHTp)^^ 

0(t) = 0(Tp) + (t - Tp) - v)> 

^(t) = ^(Tp) + (t - Tp) (^i^3(Tp) - up(Tp)) , 

where Tp is a time at which an attitude estimate is desired, v is the 
orbital true anomaly rate, and (jq's are estimated angular velocity 
components of the spacecraft system with respect to an inertial system. 
These components are estimated from the measured time history of the 
flywheel momenta, and by assuming the spacecraft is "on the average" 
local-vertical stabilized and such that the external torque is a slowly 
varying function of time. 

The transits closest to Tp which yield a strong geometry (not ill- 
conditioned) are then chosen to be used in the least squared solution. 

The resulting linear equations in cp(Tp) , 0 (Tp) , and t weighted, 

with greater weights (generally) for transits near Tp. The computer 
running time to provide one attitude estimate is approximately 0.1 seconds 
on the GDC 6400. 
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Error Analysis 

The greatest error source does not originate in the instrumentation, 
but in the attitude model. However, the estimated attitude is quite good 
even in the presence of a poorly controlled spacecraft and a poor stellar 
field. The active system has a greater accuracy than that of the passive 
system. One -sigma values of the absolute value of the error in the active 
or passive and single-head or two-head systems are: 



Two-Heads 

Single-Head 


Yaw 

Pitch 

Roll 

Yaw 

Pitch 

Roll 

Active (30 sec) 

.010° 

.014° 

.022° 

.047° 

.074° 

. 118° 

Passive 

.047° 

.027° 

.038° 

.121° 

.254° 

.367° 


These values were obtained from the poorest stellar field over all orbital 
Q's and a rather poorly controlled spacecraft attitude. 
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APPENDIX A 

ANGULAR MOMENTUM OF A BODY COMPOSED OF 
CONNECTED RIGID SECTIONS 


Assume the body is composed of two rigid sections such that these 
sections have at least one point in common (see Figure lA) . Denote the 

set of points in one rigid section as s and the set of points in the 

other as m. Assume one point of s is fixed with respect to an inertial 
system, and let; 

be an inertial system, 

be a system fixed in s at its center of mass, 

c be a vector from the center of mass of s to a common point of 

s and m, 

S be a system fixed in m at the terminus of c. 
m 



Figure lA: The Connected Rigid Bodies s and m. 
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Now, 


H = 


E R, 
T 


X m, 
1 


dt 


= E 
T 


X m 


.(■ 


d R. 
s 1 

dt 


0) 


Is 


X R 


i) ’ 


where 


T = s + m , 

H angular momentum of the body, 

R. position vector of the i^^ point in the body, the origin of R^ 
is the center of mass of s , 

f" 

is the mass of the point , 


d.R . - 

"TT = x.i. + y.j. + z.k. 
dt j j j J 


(with R = x.i. + y.j. + z.k.), 
J J J J J 


is the angular velocity of with respect to S^. 


Thus , 


H = 


S R. X 
m 


m . 
1 


d R. 
s 1 

dt 


S m.R. 

T 1 1 


Ss * 


V. 


s ince 


d R. _ 

= 0 for R.C s. 

at 1 


Also, 


Z m^R^ X X R^ = #(T, S^) , 
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where 


f(T, S ) is the inertia tensor of the total body computed with 
respect to S^. 


For CZ let = c + r^. So, 


d R. 


d c dr. 


VJ , /U V- 1. . 

T— 1 SI. / I \ /Si r I 

E R. X m. , , =E (c + r.) xm. (“J^r + "37“ ^ 

X 1 dt 1 1 \dt dt 

m m 


o) X r . 
sm 1 


= Z (c + rp X m^(u)^^x r.) 
m 


= E c X m (m x r ) + $(m, S ) m 

1 sm 1 m sm 

m 


Hence, 


H = 


HT. sp + #(m, S^) c« 3 ^ + c x(c«^^ x r J . 


where 


vector from the terminus of c (origin of S ) to the center 
of mass of m, 


M mass of m. 
m 


If several systems, m^, are present, then 


H = f (T, S ) 0 )^ + s($(m., S .) 

s Is A 1 mi 


0) .H-M.c.x(o) .xr..) 

smi ml 1 smi ®i 


If each of the systems are connected to s at their center of mass, then 
r^ . = 0. So, in this case, 

°i 


H = $(T, sp u)i 3 + Z l(m., 

1 
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APPENDIX B 

INERTIA TENSOR OF THE SOLAR ARRAY 


Figure IB represents the approximation of the dynamic configuration 
of the spacecraft which will be used to compute its inertia tensor. The 
solar array is represented by two rectangular thin plates which are both 
tipped out of the plane defined by and the rotation axis of the 
array. This tip angle, a, is used to position the panels more nearly 
in the sun’s direction. The angle, a, is a constant, but g increases 
by 360^ in an orbital period. 


To aid in the computation of the solar array inertia tensor, introduce 
a system, S , with origin at the center of mass of the right-hand panel 

P ^ A A 


as pictured 
to its edges. Then, 


Let ip be normal to the panel, and let and be parallel 



is the inertia tensor of the right-hand panel, computed with respect to 

S . Here, 

P 

total solar array mass, 
panel width, 
panel length. 

Recall that is a principal direction of the spacecraft without 
the solar array. This direction is nominally that of the pitch axis, 
which, in turn, is nominally that of the axis of rotation of the array. 
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Now the orientation of with respect to is defined by the 
angles p, p, and o as follows: 


s 


k" 


rotation -p about 
rotation p about 
rotation -a about k^' 




i" 
J? , 


= k 


Hence, 




•sin CT 
cos a 
0 



1 

0 




= S 



■ ^ • 


As before, let be parallel to S^, but with origin at the center of 

mass of the array. Then, a vector from the origin of to the origin of 
is (Figure IB) 


I j. 




— + I sin p ) iy 


— I cos pjj^ - T "23 / 


W =f 


2 2 
h2 + h3 

■'^1^2 

-^1^3 


-h^h2 

2 2 
h^ + h3 

■’^ 2^3 


-h^h3 

-h2h2 
2 2 
+ ^2 


+ S' $(R, S ) S 
P 
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-f-1 

= - (— S2j^ + -t sin p) , 
-t-l 

h2 = -(— S22 + I cos p), 


h is 

3 2 23 ■ 


The above expressions yield the inertia tensor for the right-hand 
panel. For the left-hand panel, the function is obtained by evaluation 
at -I, The same matrix, S, is assumed for both panels (that is, 

the panels are assumed to remain parallel). Hence, 


h " W W = $(m^, S^) 


®21 

■®21 ®22 

■"21 

®23 


1 

"®22 

®23 



1 - ; 

s2 

23 

2 

- "31 

■®31 ®32 

"®31 

®33 


-4 

"®32 

®33 


1 - s: 


cos pit cos p + 


— cos 0sina + -csinp cos p - sin p 

sin p(^, sin p + -^2^321^ “ T"^23 P 


2 23 


j( + cos a 
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If it is assumed that p = 0, then the above expression can be 
greatly simplified. This simplification has been used. Hence, 


Cf + cos 23 -q^ cos p -q^ sin 2p 


$ = 
s 



where 


"1 = 


^1 12 212 

— (1 - — sin a) + -t, + I 2 + I cos a 


^s ’ 


1..2 


1 .2. 


3(^1 4 V ^s’ 


1 / .2 . 2 ,1 2 . „ 

" 6^"^1 4 ^ 2 ^ ^ s ’ 


'f'l c + -| ^) M^. 
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THE SKEW- SYMMETRIC MATRIX AA* AND THE ANGULAR VELOCITY 


In the derivation of the equations of motion, it was expedient to 
trace the following path: 

(1) Write the angular momentum of the body with components 
in the body- fixed system, S^. 

(2) Transform this vector into one with components in the 
inertial system, Sj. 

(3) Set the time derivative of the transformed vector equal to 
the external torque (components in S^) . 

(4) Transform both members of the equality obtained in step 3 
into vectors with components in the body- fixed system. 

This process introduced a matrix of the form AA*, where A is a 3 x 3 
orthonormal matrix. Since the form of this matrix may not be well known, 
a short discussion concerning its relationship to angular velocity is 
now given. 


Theorem 1: 


For any position R, there exists a unique q, independent of R, such 


that 


d.R d.R 
1 _ 1 


dt dt 


+ q x R , 


where 


- 7 — = x, i, + y. j.+zk.. 
dt 1 1 *^1 L 


Here x., y. , and z. are the component of R in S., 

1 I X 1 
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Proof of Existence 
Let 



where A may be a non-constant function of time. 



Now, the matrix AA' is skew-symmetric for 

AA' = I, the 3x3 identity, 

AA' + AA' = 0 , 

AA' + (A'A) ' =0, 


Hence , 
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which proves the skew- symmetry. Hence, 



Expressing this result in vector notation completes the existence. 
Proof of Uniqueness 

Suppose there exists i- q2> then (q^ - q2) x R = 0 for all R; 
but this is impossible. This completes the proof of this classical 
theorem. 


Notation: 


q is written q = > 

with respect to S^. 


and is called the angular velocity of 


S. 

J 


Corollary: 
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where 


(JO . . (Di i . + 0)o j 4 ^ 4 • 


ij ^1 

The proof is contained in the proof of Theorem 1. 


The corollary gives a constructive method of computing the components 

of u) . . 5 but in practice this method is cumbersome. A second method will 
iJ 

now be given (the standard method) which is a result of two simple 
theorems . 


Theorem 2: 


(JD .1 - (JO . . + O) -1 . 

ik LJ jk 


Proof 


d.R d.R 

L 


So, 


dt 

dt 

d.R 

d, R 
k 

dt 

■ IT 

d.R 

L 

1 

dt 

d t 


+ (JD . . X R , 


+ (JO X R . 
jk 




Hence , 


U)ik - 0)ij + («jk 


105 




APPENDIX C 


Theorem 3 : 


If 




sin 9 
cos 9 
0 



then 


u) . . — 9 k . -9k. 
ij 1 J 


Proof 



-cos 9 
-sin 9 
0 


\ 





so the corollary implies the stated result. 


Example 


^1 ^2 notation cp about 


^2 ^3 ® about ^2 " ^3' 


What is 0)^2' 


Solution 1 


^13 ^12 ^23 


= cp + 0 ^3 ^ ® ^3 cp(sin 0 j 



+ cos 0 k^) , 


106 




APPENDIX C 


Solution 2 



sin cp 

cos 0 cos cp 
-sin 0 cos cp 




cos cp 

-cos 0 sin cp 
sin 0 sin cp 


+ 0 I sin 0 sin cp -sin 0 cos cp 

cos 0 sin cp -cos 0 cos cp 



So 



0 

0 

1 


0) 


13 


= 0 i^ + cp sin 0 + 9 cos 0 • 
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THE EXTERNAL TORQUE 
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APPENDIX D 
THE EXTERNAL TORQUE 

The formulas which are used to compute the external torques will 
now be stated without derivation. Three sources of torques are 
assumed; namely, solar pressure, gravity gradient, and magnetic. 

Solar Pressure 


T^ (solar pressure) = 2,0 x 10 


-5 



lb force - ft 


if 0 < V < 112.5 


or 247.5° < v < 472.5°, 


= 0 elsewhere. 


In deriving this expression, it was assumed the sun is at the 
equator, and sunlight only impinges upon the ^olar panels. The 
following numerical values were assumed: 


-8 2 

P = 9.4 X 10 lb force/ft (solar pressure constant), 
6o = . 1 (panel albedo), 

H = 4 ft, 

= 2 ft, 

.(,2 = 8 ft (see Figure IB) • 

Gravity Gradient 


(gravity gradient) 



cos 2p)cp - q 2 0 cos p 
2q^^ cos 2p)0 - q 2 cp cos p 
q 2 cp cos p - q 2 9 sin p 


q 2 sin p 
q^^ sin 2p 
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where 

-M 2 = = 3,000 pole-cm, 

-8 

K = 1.477 X 10 lb force - ft/pole-cm. 

The expression was obtained from an unpublished G. E. document. 
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APPENDIX E 

DETERMINATION OF APPROXIMATE INGRESS, EGRESS, AND TRANSIT TIME 

As an aid to the star identification, it is desirable to compute 
the approximate time any star enters the field of view (ingress) , leaves 
the field of view (egress), and transits the slit. These approximate 
times may be obtained by assuming the spacecraft has its nominal 
attitude, that is, roll, pitch, and yaw are zero. 

In this study, each star is assumed identified, so the identification 
problem is not considered. The approximate ingress, egress, and transit 
times are needed since initial estimates are required to obtain the exact 
transit times (necessary for a simulation since measurements are not made) 
via a recursive method. 

Approximate Transit Time 

If\j/=0=cp = O, then from Equation (16) 

0 = B^(y)S(a,3)B(-v)h 

at the instant the star transits the slit (the subscripts i and j have 
been dropped) . This equation is simply one equation In the unknown 
transit time. The true anomaly, v, and the slit orientation angle, 
are given functions of time. So, 

/ cos (v-p)h^ + sin (v-p)h^ \ 

0 = (cos Y, 0, sin y)S(q,0) ( ^2 ] 

\-sin(v-p)h^ cos (v-p)h^ J 

or 

0 = a cos (v-p) + b sin (v-6) + c, 


where 
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a = cos Y - sin Y <^°s a h^, 
b = sin Y cos a + cos y 
c = sin Y ^ ^3 ^ ^^3* 


Hence , 


-bvh2 + ^ V 1 

sin (v-3) = 5 5 , 

(1 - hp (1 - v^) 

-avh- + b Vl-v^-h? 

cos (v-p) 5 5 . (2E) 

(1 - hp (1 - vp 

A problem remains as to which of the above signs are to be used in 
Equation (2E) . This choice may be resolved by determining cos T| and 
sin T| from Equation (15). Hence, 


/ cos (v-3)h^ + sin 
cos Tj = B^SCqjO) j Y 12 

\ -sin (v-p)hj^ + cos 


(v-p)h2 

(v-p)h2 


or 


(l-vp 


cos T| = h 2 cos a + sin a cos y 



Similarly, 


(l-vp 


sin T| = -h^ 


cos Y sin a + cos a 



Now all transits must be such that 0 < T| < 90^; also, the cant 

angle, a, will be chosen such that 0 < a < 90^. So, if cos y ^ Oj the 

lower sign must be chosen. If cos y < 0, the upper sign must be chosen 

if a > r^=T! (one-half the field of view). 

2 max 
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Hence, if y is a given constant (passive system), the approximate 
time of transit may be obtained as follows: 

-h w ± a 

t = to + T tan“ 

V - r 

-a V h 2 + b V 1 

where the upper sign is chosen if cos y > and the lower sign if cos y < 0 
and a > r 2 * The quadrant is chosen by the sign of the numerator and 
denominator. A circular orbit is assumed (that is, v = constant). 

If y is not given (active system) , then Equation (3E) is not a 
solution for the transit time. In this case, we assume 

y = yo + yt, 

where y^ and y are given constants and y > > v. Equation (IE) may be 
written 

0 = A cos y + B sin y, 

where 

A = cos (v-p) + h^ sin (v-p) , 

B = cos a ("h^ cos (v-p) + h^ sin (v-p)) + sin a h^ . 

Hence , 

+ A 

sin y = 

+ B 

cos y = 
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But , 


sin n = + vA^ + , 


therefore , 


tan Y = — 


Thus, the transit time may be found from the recursive relationship 


+ 4 tan-1 


A(ti> 

-B(t.) ’ 


which converges for y > > v. 


Approximate Ingress and Egress 


At ingress or egress, the angle between the star and is r^j the 
half field of view angle. Hence, 


cos ^2 “ sin a) 


cos (v-R)hj^ + sin (v“3)h^ 


-sin (v-p)h^ + cos 


So, at an ingress or egress, 


sin (v“3) = 


-hj^ c ± ^3 V (p-h2> (h2~q) 
2 

sin a (l-h 2 ) 


cos (v-3) = 


-h^ c + hj^ V (p-h 2 > (h 2 -q) 
2 

sin a(l - h 2 ) 
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where 


c = cos a h 2 - cos r 2 , 

p = cos (a + r 2 > , 

q = cos (a - r 2 > . 

o ^ ^ 

Note from Figure 2 that for 0 < CT < 90 , • h < 0 at ingress, 

/V ^ 

and ig • h > 0 at egress. But, 

± V (p-h2> (h2-q) 

^9 ^ sin a 

So, choose lower sign for ingress, upper for egress. 
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